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Abstract. We introduce a new map from polynomials orthogonal on the 
unit circle to polynomials orthogonal on the real axis. This map is closely 
related with the theory of CMV matrices. It contains an arbitrary parameter 
A which leads to a linear operator pencil. We show that the little and big 
-1 Jacobi polynomials are naturally obtained under this map from the Jacobi 
polynomials on the unit circle. 



1. Introduction 

In this paper we propose a new map from polynomials orthogonal on the unit 
circle (OPUC) to polynomial orthogonal on the real line. This map can be consid- 
ered as a one-parameter generalization of the Delsarte-Genin [19] map from OPUC 
to symmetric orthogonal polynomials on an interval. We derive this map from ele- 
mentary properties of the CMV matrices. The CMV matrices are known to play a 
crucial role in the theory of OPUC [31], [35], [11]. 

This new map contains an additional arbitrary parameter A which leads to infin- 
itely many families of polynomials P n (x; A) orthogonal on the real line. In general, 
for arbitrary real values of A these polynomials are orthogonal on a union of two 
distinct intervals of the real line. 

Starting from given explicit sets of OPUC, one thus constructs families of poly- 
nomials orthogonal on the real line P n (x; A). The main result of this paper consists 
in the following. Starting from the Jacobi polynomials on the unit circle, we con- 
struct a one-parameter family of polynomials P n (x; X) orthogonal on the union of 
two intervals of the real line. 

We identify these polynomials with the big -1 Jacobi polynomials discovered 
recently in [38] . For a particular value A = 1 of the parameter we obtain the little 
-1 Jacobi polynomials found in [37]. 

The underlying spectral problem has the following form 

(1.1) (L + AM — xl)q(x, z) = 0, 

where L and M are tridiagonal matrices of a special form (L and M have zeros on 
the super- and sub-diagonals). 

If we fix A, the problem (1.1) is a classical eigenvalue problem for the tridiagonal 
matrix K(X) = L + AM. Thus, it gives rise to a system of orthogonal polynomials, 
Pade approximants and moment problems. 
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If wc fix x, the problem (1.1) is a generalized eigenvalue problem for the two 
matrices M and L — xl . It generates a system of functions rational in x that are re- 
lated to the theory of bi-orthogonal rational functions [43] . Moreover, it has recently 
been shown that such generalized eigenvalue problems appear in rational interpola- 
tion [6], [18], Nevanlinna-Pick problems [17], [18], integrable systems [33], [36], and 
mechanics [24, Section 2.4] (see also [20, 21]). 

A priori, the general theory associated to the problem (1.1) is rather involved. 
That is why we make some assumptions on the matrices L and M. Namely, we 
assume that L 2 = I and M 2 = I. Thus, the main idea is to identify systems of 
orthogonal polynomials that have a hidden biorthogonality associated to a param- 
eter. 

Note that the special case x — leads to the Schur linear pencil L + AM con- 
sidered by Watkins [35]. 

Finally it should be mentioned that particular cases of (1.1) surprisingly arise in 
the theory of martingale polynomials [9], [30], two-state free Brownian motions [1], 
and Gaussian processes [27]. Let us also signal references [22] and [16] where sim- 
ilar pencils appear in the investigation of exactly solvable kinetic models and the 
corresponding combinatorial problems. 

The paper is organized as follows. 

In Section 2, we recall elementary facts concerning OPUC and CMV matrices. 

In Section 3, we introduce the Jacobi (3-diagonal) matrix J starting from two 
unitary elementary matrices L and M. A new family Q n (x) of orthogonal polyno- 
mials on the real line is associated with the matrix J. 

In Section 4, basic results on the Delsarte-Genin map from symmetric orthogonal 
polynomials S n (x) to OPUC are presented. 

In Section 5, we introduce the companion Delsarte-Genin polynomials P n (x) as 
Christoffcl transforms of the polynomials S n (x). We show that the polynomials 
Q n (x) coincide with P n (x). 

In Section 6, we consider the Jacobi OPUC and show that the corresponding 
polynomials Q n (x) coincide with the little -1 Jacobi polynomials. 

In Section 7, a one parameter generalization of the previous scheme - namely the 
Schur-Delsarte-Genin (SDG) map - is proposed. Starting from the same unitary 
matrices L and M we construct a 3-diagonal matrix K(\) = L + AM that depends 
on a parameter A. The matrix J corresponds to A = 1. The matrix K(X) generates 
a new family of orthogonal polynomials Q n {x;X) depending on the parameter A. 
A linear operator pencil is associated with the eigenvalue problem for the matrix 
K(X). 

In Section 8, we consider the special case oo = Oi = ■ ■ ■ = that corresponds to 
the OPUC $> n (z) — z n . The polynomials Q n {x\ A) are then shown to be orthogonal 
on the union of two intervals of the real axis. 

In order to recover the orthogonality measure for the general case, we need 
a special map from symmetric to non-symmetric orthogonal polynomials that has 
been proposed by Chihara. This leads to a nontrivial linear operator pencil J±+\J2- 
This technique is described in Sections 8 and 9. 

In Section 10 we relate the polynomials Q n (x;\) to special polynomials P n (x) 
independent of A and having a known orthogonality measure. 
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Finally in Section 11, we apply this technique to Jacobi OPUC. The correspond- 
ing polynomials Q n (x\ A) are shown to coincide with the big -1 Jacobi polynomials 
which are orthogonal on the union of two intervals on the real axis. 

2. Polynomials orthogonal on the unit circle and CMV matrices 

Monic polynomials $„(z) = z n + 0(z n ^ 1 ) on the unit circle are defined through 
the recurrence relation 

(2.1) <f> n+1 {z) = z$ n (z)-a n $* n (z), 
with initial condition &q(z) = 1 where 

$;(z) = z n $JT/T). 

The recurrence coefficients a n = — $„ + i(0) are called reflection parameters (some- 
times also referred as the Schur, Geronimus, Verblunsky, ... parameters). It is well 
known that the reflection parameters are characterized by the property \a n \ < 1 for 
n = 0,1,2,.... 

In what follows we will consider only the case where a n are real. 

For convenience, we set 

(2.2) a_i = -l. 
The dual recurrence relation is 

(2.3) $* n+1 {z) = $>* n {z)-a n z$ n {z). 
It is standard to introduce the parameters 

r n = V 1 - a l 

(the arithmetic meaning of the square root is assumed). 

Following [35] and [11], let us introduce the block-diagonal matrix 



(2.4) 



L = 



ro 



V 

and the block-diagonal matrix 



a 2 



T2 
-0,1 



(Z4 



— (24 



/I 



(2.5) 



M 



ai 



n 
-oi 



a 3 
r 3 



r3 

-a 3 



a 5 



-a 5 
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Both L and M are block-diagonal unitary matrices. Define also the 5-diagonal 
unitary matrix U [35], [11] 

(2.6) U = LM 

(which is called the "Zigzag matrix" in [35] and the "CMV matrix" in [31]). 
In [28] it was proposed to consider the Hermitian matrix 

(2.7) H = U + Ul 

The authors of [28] connected this matrix to a direct sum of two Jacobi matrices 
and then arrived at the famous Szego relations between polynomials orthogonal 
on the unit circle and on the real line. In the next section we propose a slightly 
different approach to the matrix H which connects the polynomials on the unit 
circle with a more simple class of polynomials orthogonal on the interval. These 
polynomials are related to symmetric orthogonal polynomials on the same interval 
by a simple Christoffel transform. 

3. From CMV matrices to polynomials orthogonal on the real line 
On the natural basis e„, n = 0, 1, 2, . . . the matrix H is 5-diagonal and Hermitian 
He n = r n r n+ ie n+2 + r n (a n+ i - a n _i)e n+ i - 2a„a„_ie„ + 

(3.1) r n -i(a n - a n _ 2 )e„_i + 

Introduce the 3-diagonal matrix J which is defined by 

(3.2) Je„ = r n e n+ i + (a n - a„_i)e„ + r n _ie„_i. 
It is easily verified that 

(3.3) H = J 2 - 2. 

Thus among the eigenvectors of the matrix H there are the eigenvectors of the 
matrix J. Moreover, it is seen that 

(3.4) J = L + M. 

Note that the matrices L and M are both unitary and satisfy the obvious property 
L 2 = M 2 = 1. Hence 

(3.5) J 2 = (L + M) 2 = L 2 + M 2 + LM + ML = 2 + U + U*=2 + H, 

thus recovering relation (3.3). 

The matrix J is a symmetric Jacobi matrix. It generates orthonormal polyno- 
mials Q n {x) satisfying the recurrence relation 

(3.6) r n+1 Q n+1 (x) + (a n - a„_i)Q„(x) + r n Q n -i(x) = xQ n (x). 
The corresponding monic orthogonal polynomials 

(3.7) Q n (x)= ® n{x) =x n + 0(x n ~ 1 ) 

T\T2 ■ ■ ■ r n 

satisfy the recurrence relation 

(3.8) Qn+i{x) + b n Q n (x) + u n Q n _i{x) = xQ n (x) 
with initial conditions 

Qo = 1, Qi(x) = x-bo, 
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where the recurrence coefficients are 

(3.9) u» = r£_i = 1 - al-i, b n = a n -a n -i. 
Note that 

u = 0, b Q = a + 1 

due to the assumption a_! = —1. 

The relation between z and x is given by 

(3.10) x = z l / 2 + z- 1 ' 2 . 

Here it is assumed that one branch of the two-valued analytic function z 1 ! 2 is 
chosen. In what follows we shall take z 1 / 2 = r 1 / 2 e^/ 2 when z = re 1 '*', < <j> < 2ir. 
Suppose that a point z is moving counterclockwise on the unit circle \z\ = 1 starting 
from the point z = 1 and returning to the same point after a full rotation. Then 
the corresponding image point x — z 1 / 2 + z~ x l 2 is moving in one direction on the 
interval [—2, 2] starting from the point 2 and ending at the point —2. Thus the 
mapping (3.10) is a one-to-one correspondence between the punctured unit circle 
(i.e. z / 1) and the interval [—2,2]. The only exception is the point z = 1: it 
corresponds to the two endpoints ±2 of the interval. 

Generic polynomials <& n {z) under condition \a n \ < 1 (with complex coefficients 
a„) arc orthogonal on the unit circle 

(3.11) / § n {e ie )$> m {e- w )dv(0) = 0, n + m, 
Jo 

where v{0) is a nondecreasing function in the interval [0, 27r]. 

For the case of real parameters o„ we will assume that dv(6) — p(9)d0, where 
p{9) is a positive weight function: 

(3.12) / $ n (e w )<S> m (e- M )p(e)de = 7 n^m. 
Jo 

The function p(8) satisfy the symmetric property 

(3.13) p(2tt -6)= p{6) 

which is valid for real coefficients a n verifying the condition — 1 < a n < 1 [31]. 

To end this section, note that Q n (x) are polynomials orthogonal with respect to 
a measure. In fact, this measure can be found by means of the map proposed by 
Delsarte and Genin. 

4. The Delsarte-Genin map to symmetric polynomials on the interval 
Consider the symmetric polynomials S n (x) satisfying the recurrence relation 

(4.1) S n+1 (x) +v n S n - 1 (x) = xS n (x), 
where 

(4.2) u„ = (l + a„_i)(l-a n _ 2 ) 

Delsarte and Genin showed [19] that the polynomials S n (x) are related to the 
polynomials $„(z) as follows (see also [40]) 
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The reciprocal map is 

z"/ 2 (z 1 /2 5n+l(a . ) _ CTn 5 n ( a . )) 



(4-4) 



z-l 

z (n+l)/2 {z l/2 Sn+i{x) _ anZ l/2 Sn{x)) 



1 - z 
where 

The reflection parameters a n are expressed as 

(4.5) fln-l = 1 - cr„. 

The symmetric orthogonal polynomials S n (x) are orthogonal on the interval [—2,2]: 

(4.6) J S n (x)S m (x)w(x)dx = 0, n ^ to, 

where the weight function u)(x) is symmetric w(—x) = w(x) and is related to the 
weight function p{9) in the orthogonality relation (3.12) as follows [19, 40] 

(4.7) "M = **m> * = *™m- 

Note that sin(6*/2) > for < 9 < 2n and that the symmetry property (3.13) is 
equivalent to the symmetry property w(— x) = w(x). 

There is an obvious "scaling" generalization of the Delsarte-Genin mapping [40]. 
For an arbitrary positive parameter g define the map 

(4.8) ^fert- ^ 1 ( * >W + * iW) . 

1 - a„_i 

where the relation between x and z is 

(4.9) x = g(z 1 / 2 + z- 1 / 2 ). 

It is easily seen that the polynomials S n (x; g) are symmetric orthogonal polynomials 
satisfying the recurrence relation 

(4.10) S„(x;g) + g 2 v n S n -i(x;g) = xS n (x) 

with v n also given by (4.2). The polynomials S n (x) introduced in (4.3) correspond 
to the case g = 1. From (4.10) it is clear that 

(4.11) S n (x;g) =g n S n (x/g) 

i.e. the polynomials S n (x;g) are obtained from the polynomials S n (x) by a scaling 
transformation of the argument x. In particular, if the polynomials S n (x) were 
orthogonal on the interval [—2,2], the polynomials S n (x;g) are orthogonal on the 
interval [2g,—2g]. 

The reciprocal transformation is given by the formulas 

(4.12) <p n (z)=g — , 

z — l 

where 

S n+1 (2g;g) 
ani9) = ^J2gW = 9 ^ 
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and 

5. The companion Delsarte-Genin map and orthogonality measures 
for sdg-transformed polynomials 

In this section we find the orthogonality measure for Q n (x) by performing a 
Christoffcl transformation of the polynomials from the previous section. The result- 
ing map will be called the Schur-Delsarte-Genin map (or, shortly, the SDG-map). 

First, consider the polynomials P n (x) which are obtained from S n (x) by the 
Christoffel transform 

(5.1) P n { x ) = Sn+1 -^f n{x \ 

X ~ p A 

where 

(5.2) A n = = -o n = a„_i - 1 

(the last formula follows from the property S n (—x) = (— l) n S n (x) for the symmetric 
orthogonal polynomials). 

The polynomials P n {x) satisfy the recurrence relation 

(5.3) P n +i{x) + b n P n (x) + u n P n -i(x) = P n {x), 
where (see, e.g. [42]) the new recurrence coefficients are 

A 

(5.4) u n = v n , b n = A n+1 -A n . 

From (5.4) and (4.5), we immediately see that the recurrence coefficients b n ,u n 
coincide with (3.9): 

u n — 1 — a n -ij b n — a n — a n —\. 

This means that the polynomials P n (x) and Q n {x) coincide. Thus the (non- 
symmetric) polynomials Q n (x) are obtained from the symmetric polynomials S n (x) 
by the Christoffel transform (5.1). This allows expressing the polynomials Q n (x) 
in terms of the circle polynomials & n (z): 

(5.5) w^- ^-y' 1 

The reciprocal transformation from the polynomials S n (x) to the polynomials 
Q n (x) is given by the Geronimus transform 

(5.6) S n (x) = Q n {x) - B n Q n _ 1 (x), 
where 

(5.7) B n = 2- A n = l + a n - 1 . 

The polynomials Q n {x) are orthogonal on the interval [—2, 2]: 

(5.8) j Q n (x)Q m (x)w(x)(x + 2)dx = 0, n 7^ to, 

where the function w(x) is given by (4.7). Formula (5.8) follows from the definition 
of the Christoffel transform [34] . 
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Consider the special value Q n (2). From (5.5) we have (x = 2 corresponds to 
z = l) 

(5.9) Q n (2) = $ n (l) = (1 - ao)(l - ai) ... (1 - o„_i). 

Using these formulas, one can easily obtain an expression for the polynomials $ n (z) 
in terms of the polynomials Q n {x): 

f r lnl ft / x „n/2 (^' 1/2 - an)Qn(x) - Q n +l{x) 

(5.10) $„(z) = z/ _ . 

Note that apart from the polynomials Q n (x) one can introduce the adjacent poly- 
nomials Q^n\x) which are obtained from S n (x) by the Christoffel transformation 
which is adjacent with respect to (5.1): 

(5.11) QV{ X )= Sn+i ~ x A :l Sn{x \ A( n ) =^ l =-^= 1 -^- 

The monic polynomials Qn ■* (x) satisfy the recurrence relation 

(5.12) Q ( -^(x) + b^Qi~\x) + u^Q { -\{x) = xQ^\x), 
where 

(5.13) b ( - } = a n _i - a n = -b n , = 1 - a 2 n _ x = u n . 

Comparing (5.12) with (3.8), we see that the polynomials Q < k \x) have the same 
recurrence coefficients u n as the polynomials Q n {x) but with coefficient b n that 
have the opposite sign. This means that these polynomials almost coincide: 

(5.14) Qi~\x) = {-\TQ n {-x). 

The polynomials Q\i\x) are orthogonal on the interval [—2, 2] 

r 2 

(5.15) J Q ( - ) (x)Q^(x)w(x)(x-2)dx = Q, n^m. 

We see that the polynomials Q n {x) are more convenient to construct the polyno- 
mials $„(z) than the polynomials S n {x). Indeed, the reflection parameters a n are 
immediately determined, up to a sign, from the recurrence coefficient u n = 1 — a 2 l _ 1 : 

(5.16) a n = ±y/l - u n+ i. 

Choosing the sign in (5.16) selects the polynomials Q n {x) or the polynomials 

Q ( n\x). 

We can consider a third set of orthogonal polynomials T n (x) which is obtained 
from the polynomials Q n {x) by the Christoffel transform at x = 2: 

T ( \ Qn+l(X) 9 ^ lQn{x) 9^±M ZSLZ a r,)Qn{x) 

(5.17) l n [x) - — - — . 

Simple calculations show that the polynomials T n (x) are symmetric orthogonal 
polynomials satisfying the recurrence relation 

(5.18) T n+1 (x) + w n T n -i(x) = xT n (x), T = 1,T 1 = x, 
where 

(5.19) w n = (1 + a„_i)(l - a n ). 
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Note that the polynomials T n {x) can be obtained from the polynomials S n (x) by 
a "symmetric Christoffel transform" (see, e.g. [34], [33]) 

T = S n+2 (x) - ?$±ffis n (x) = Sn+2 (x) - (1 - a n )(l - an-^S^x) 

(5.20) " ^-4 x 2 -4 

_ xS n+ i(x) - 2(1 - a Tl _i)5 n (x) 
x 2 - 4 

Using formula (4.3) and the recurrence relations for the polynomials <&„(z) we can 
express the polynomials T n {x) in terms of the polynomials & n (z): 

(5.21) T n (x) = *- /2 (**n(*) -*;(*» 

z — 1 

The polynomials T„(a;) are orthogonal on the interval [—2, 2]: 

(5.22) J T n (x)T m (x)(4 - x 2 )w{x)dx = 0, n^m. 

Thus with the polynomials Q n (z) orthogonal on the unit circle with real reflection 
parameters \a n \ < 2 one can relate 3 types of polynomials S n (x) , Q n (x) , T n (x) , 
orthogonal on the interval [—2, 2] of the real line. Two of these polynomials, S n (x) 
and T n (x), are symmetric, whereas the polynomials Q n {x) are non-symmetric. 

Note that the polynomials S n (x) and T n (x) coincide with the so-called "com- 
panion polynomials" considered in [8]. Formulas (4.4) express the polynomials 
$„(z) in terms of the polynomials S n (x). Note also that the relations between the 
symmetric polynomials T n (x) and the non-symmetric polynomials Q n {x) (without 
the identification of the connection to OPUCs) were studied by L.M. Chihara and 
T.S. Chihara in [15]. 



6. From the Jacobi OPUC to the little -1 Jacobi polynomials 

The unit circle analogues of the Jacobi polynomials were proposed by Szego [34] 
and studied by Golinskii [25] and Badkov [4], [5]. These Jacobi OPUC $ n (e ie ) arc 
orthogonal on the unit circle with respect to the weight function 

(6.1) p{6) = (l-cos0)« +1/2 (l + cos( 



^+1/2 



where £ [0, 2ir}. It is assumed that £ > —1, i] > —1. 
Their parameters a n have the following expression [4] 

f , T iT^ i n if n is even 



(6.2) 



- if n is odd 



The condition a_i = —1 obviously holds. 

Note that the Jacobi OPUC are connected with the generalized Gegenbauer 
(ultraspherical) polynomials through the DC map. 

Indeed, the generalized Gegenbauer polynomials [14], [7] are symmetric orthog- 
onal polynomials satisfying (4.1) with 



(6.3) 



n(n+2Q 
(n-K+J7)(n-K+rH-l) 



if n is even 



if n is dd 
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where £ > — 1, r/ > — 1. These polynomials are orthogonal on the interval [—2, 2]: 

(6.4) J S n (x)S m (x)(4-x 2 )< i \x\ 2r >+ 1 dx = 0, n^m. 

Using the Delsarte-Genin map (4.4) and formula (4.2), we can easily verify that 
the Jacobi OPUC $„(e* e ) correspond to the generalized Gegenbauer polynomials 

S n (x). 

Consider now the companion polynomials Q n (x) corresponding to the parame- 
ters (6.2). They are orthogonal on the interval [—2, 2] with respect to the function 

w(x) = (x + 2)(4-x 2 )«|x| 2 " +1 . 

The recurrence coefficients b n , u n for these polynomials are explicitly obtained from 
formulas (3.9), where the coefficients a n are given by (6.2). 
Choose the parametrization 

, a - 1 P-l 



2 ' ' 2 

A simple analysis shows that the polynomials Q n (x) coincide (up to an obvi- 
ous affine transformation of the argument) with the little -f Jacobi polynomials 
Pn _1 \x; a, f$) introduced in [37]. The little -1 Jacobi polynomials possess a re- 
markable property: they are polynomial eigenfunctions of a Dunkl-type differential 
operator of the first order. Hence these polynomials provide a new nontrivial ex- 
ample of "classical" orthogonal polynomials on the real line (for details sec [37]. 

Note that similar polynomials were considered in [15], where the recurrence coef- 
ficients b n , u n were calculated explicitly. These polynomials belong to a special class 
of orthogonal polynomials and were studied by another approach in [2]. Our ap- 
proach yields explicit expressions for the recurrence coefficients and the orthogonal 
polynomials Q n (x) by a straightforward application of the correspondence formulas 
between polynomials orthogonal on the unit circle and on the interval. 

7. A ONE-PARAMETER GENERALIZATION OF THE COMPANION POLYNOMIALS 

Consider a linear pencil of matrices 

(7.1) K(X) = L + XM 

with an arbitrary parameter A. 
We have 

(7.2) K 2 {\) = L 2 + \ 2 M 2 + X(LM + ML) = 1 + A 2 + A(J 2 - 2) = 1 + A 2 + XH. 

We see that the Hermitian matrix H can be expressed in terms of K 2 (X) for every 
A. On the other hand, the matrix K(X) is a nondegenerate symmetric tri-diagonal 
matrix 

/a + A r \ 

ra -oo + Aai Ari 

Ari a 2 - Aai r 2 

r 2 -a 2 + Aa 3 Ar 3 



K(X) 



\ ■ J 

Here, nondegenerate means that all off-diagonal entries of the matrix K(X) are 
nonzero if A ^ and |a,| ^ 1. 
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Hence one can define a family of formal monic orthogonal polynomials Q n (x; A) 
depending on the argument x and an additional parameter A. These polynomials 
are uniquely defined through the 3-term recurrence relation 

(7.3) Q n +i(x; A) + b n (X)Q n (x; A) + u n (X)Q n - 1 (x; A) = xQ n (x; A), 

where 

a n — Aa n _i if n is even 
Aa„ — a„_i if n is odd 



(7.4) b n (X) = 
and 

(7.5) u n (X) = 



A 2 (l — dn-i) ^ n i s even 
1 — a^-i ^ n i s oc id 



Note that the eigenvalue problem for the orthogonal polynomials Q n (x;X) can be 
presented in algebraic form as 

(7.6) (L + XM -xl)q = 0, 

where q is a vector constructed from the (non-monic) polynomials Q n (x; A). Equa- 
tion (7.6) contains two parameters A and x and belongs to the class of the so-called 
multi-parameter eigenvalue problems [3], [32]. 

We already considered the case A = 1 that leads to the orthogonal polynomials 
Q n (x) with a positive measure on the interval [—2, 2]. For general complex values of 
the parameter A it is still unclear what the spectral properties of the corresponding 
orthogonal polynomials are because the matrix K(X) is non-Hcrmitian. Neverthe- 
less, when A is real, the polynomials Q n (x; A) are orthogonal on the real axis with 
a positive measure. Indeed, as seen from (7.5) the recurrence coefficients u n are 
positive for all n > and for all real values of A. In the special case A = — 1, 
there is a simple relation between the polynomials Q n (x; —1) and the polynomials 
(— l) n <fr n (— z) on the unit circle. 

Indeed, in the case A = —1, the recurrence coefficients of the polynomials 
Q n (x; — 1) become 

(7.7) b n (-l) = (-l) n {a n + a n -i), u n (-l) = 1 - a^. 

The coefficients (7.7) are obtained from the coefficients (5.13) by a simple trans- 
formation a n — > (— I) rl+1 a„ corresponding to the monic polynomials $„(z) = 
(— 1)™$„(— z) which are also polynomials orthogonal on the unit circle (since ob- 
tained from a simple reflection of the argument z with respect to the real axis). 
The moments c„ corresponding to the polynomials & n (z) are c„ = (— I)™c„. 
It is easily verified that 

and hence, one can introduce new symmetric polynomials T n (x) obtained from 
Q n (x; —1) by the Christoffcl transform 

* , \ _ Qn+i(x; -f ) - A n Q n (x; -I) 
ln(X) ~ x + 2 

These polynomials satisfy the recurrence relation 

f n+ i{x) + v n T n -i(x) = xf n (x) 
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with 

t5„ = (l + (-l) n a n )(l + (-l) n o„_i). 
Note that Watkins proposed [35] to consider the unitary linear pencil problem 

(7.8) (L + XM)v = 

which leads to Szcgo polynomials orthogonal on the unit circle. In this case x = 
and the parameter A belongs to the unit circle |A| = 1. 

8. The 1-periodic case 

In order to better understand the nature of the polynomials Q n (x; X) let us 
consider the simplest periodic case, that is, ao = ai = ■ ■ ■ = implying that 
ro = r\ = ■ ■ ■ = 1. The corresponding OPUC are simple monomials [31] 

(8.1) $n(z) = Z n . 

The symmetric polynomials S n (x) obtained by the DG transform (4.3) coincide 
with the Chebyshev polynomials of the first kind 

(8.2) So = l, S n (x) = 2cos(m9/2), n = l,2,..., x = 2 cos(6»/2). 
They are orthogonal on the interval [—2,2] with the weight function 

1 



(8.3) 



Wq(x) = 



Vt 



In order to recover the orthogonality measure for the polynomials Q n (x, A) we need 
some additional tools. Namely, let us introduce the following function 

(8.4) m{z,X) = {(z-K{X))- 1 e ,e ), z e C+, 

where e — (1, 0, . . . ) T . It is called the m-function (or Weyl function) of K (A) [31, 
Section 1.2]. It is a standard fact that 



(8.5) 



m(z, A) 



1 



z — A + mi(z, A) ' 
where mi(z, A) is the m-function of the matrix 



A e 



(8.6) 



/o 

A 



V 



A e 



) 



Let us also consider the following family of matrices 



(8.7) K per (X) 



1 



V 



I— 1 







A 


A 







I— 1 



Due to 
(8.8) 



3.5), we have that 



m(z, A) 



K(X) — diag(A, 0, 0, . . . ), Xe 



.(z,X) 



1 + Xm per (z, A) ' 
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Next, observe that m per satisfies the following equation 

1 



(8.9) m per (z) = - 



1 

z — 



z + X 2 m per (z, A) 
which is equivalent to the algebraic equation 

(8.10) A 2 zm^ r + (z 2 - 1 + A 2 ) m per + z = 0. 
Therefore, we have that 

(8.11) m per (z, A) = (-z 2 - A 2 + 1 + \/(z 2 - A 2 - l) 2 - 4A 2 ) , 

where the branch of the square root is chosen such that m per (z, A) ~ — - as z — > oo. 
Recall that the spectrum a(K) of the operator K is the set of all complex numbers 
( for which K — £J is not invertible. The essential spectrum a ess (K) is defined as 
<j ess {K) = <t(K) \ <Jd(K), where <Td(K) is the set of all isolated eigenvalues with 
finite multiplicity. 

The Floquet theory for semi-finite periodic self-adjoint Jacobi matrices [23] says 
that the essential spectrum of K per (X) admits the representation 

a ess (K per (X)) ={z£l:(z 2 -(l + A 2 )/A e [-2, 2])} = 
= [-A-l,-|A-l|]u[|A-l|,A + l]. 

Moreover, K per (X) has as an eigenvalue for A > 1 [23]. So, according to Weyl's 
theorem we have that a ess (K(X)) — <r ess (K per (X)) since K — K per is a self-adjoint 
one-dimensional operator. Furthermore, a{K(X)) — a ess (K(X)) because m(z, A) 
does not have isolated poles. Thus, from (8.12) we clearly see what the A-dynamics 
of the spectrum of K(X) is. Namely, the spectrum starts with two points 1 and 
-1 at A = 0. Then, the points become two different intervals moving towards each 
other while A runs from to 1. At the time A = 1, the intervals touch each other. 
After that, they are moving towards infinity in the opposite directions as A — > oo. 

Finally, substituting (8.11) into (8.8) and applying the Stieltjes-Perron inversion 
formula, we arrive to the fact that the polynomials Q n (z,X) are orthogonal with 
respect to the weight function [12] 



(8.13) w(t) -- 

In particular, for A 



\M 2 


-(i 2 -A 2 -l) 2 


-\(t 2 - 


-2A 2 i+A 2 -l) 


V 4A2 - 


-(t 2 -A 2 -l) 2 


A(t 2 - 


2A 2 t+A 2 -l) ' 


: 1 we 


get that 



te [|A-i|,A + i], 
te r— A — l, — |A — ll 



J A - f2 [2+1 

(8.14) w(t) = ^—1- = ^_±f = Wo (t)(t + 2). 

From (8.14) it is seen that the polynomials Q n {x) = Q n (x;l) coincide with the 
Chebyshev polynomials of the third kind [29] . Similarly, for A = — 1 the polynomials 
Q n (x;—1) coincide with the Chebyshev polynomials of the fourth kind. 

In fact, the behavior of the set of orthogonality is typical for the entire class of 
general matrices K{X). Moreover, the weight functions have a similar structure as 
well. However, before dealing with the general case we have to develop a special 
technique. 
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9. From the Chihara map to the generic SDG-maps 

Here, we recall some well-known results in the theory of orthogonal polynomials. 
The main goal of this section is to elaborate the machinery for the SDG-maps in 
the case of arbitrary real A, i.e. the generic SDG-maps. 

Let us start with the following simple lemma 

Lemma 9.1. Let P„(x) be monic orthogonal polynomials satisfying the recurrence 
relation 

(9.1) Pn+l(x) + {0-A n - C n )P n (x) + C n A„_lP„_l (x) = xP n ( X ), 

with the standard initial conditions 

(9.2) P = l, P 1 (x)=x-6 + A . 

Assume that the real coefficients A ni C n are such that A„_i > 0, C n > 0, n = 
1,2... and that Co = 0. Take 9 to be an arbitrary real parameter such that P n {9) ^ 

for n = 1, 2, 

Define the new polynomials 

(9.3) p n{x) = Pn +1 {x)-A n P n {x) 



x-6 

Then the monic polynomials P n (x) are orthogonal and satisfy the recurrence relation 

(9.4) P n +i(x) + (0-A n - C n+1 )P n (x) + C n A n Pn-i(x) = xP n (x). 

The inverse transformation from the polynomials P n (x) to the polynomials P n (x) 
is given by the formula 

(9.5) P n {x) = P n (x) - C n P n -i(x). 

To prove Lemma 9.1 it is sufficient to observe that relation (9.1) provides us with 
the explicit LU- factorization of the monic Jacobi matrix J — 9 corresponding to the 
polynomials P n (x) . The rest is an obvious reformulation of the theory of spectral 
(or Darboux) transformations [42], [10]. 

Note also that 

which can easily be verified directly from (9.1). 

Remark 9.2. If the polynomials P n {x) are orthogonal with respect to a weight 
function w{x) on the interval [a, [3] then the polynomials P n {x) are orthogonal with 
respect to the weight function (x — 9)w(x) on the interval \a,0\. 

The next Lemma will be useful in the identification of polynomial systems with 
known families of orthogonal polynomials. 

Lemma 9.3 ([13]). Let P n (x) and P n {x) be two systems of orthogonal polynomials 
defined as in the previous Lemma. Also, assume that the polynomials P n {x) are 
orthogonal with respect to a weight function w{x) on the finite interval [a, j3] and 
take any 9 G (— oo,a). Finally, define the following monic polynomials 

(9.7) S 2n (x) = P n (x 2 +a-c 2 ), S 2 n+i(x) = (x- X )Pn(x 2 +a-c 2 ), 
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where c is a positive number such that \J a — 8 < c and x is a real number defined 
by the relation 9 — x 2 + ct — c 2 . Then the polynomials S n (x) are orthogonal with 
respect to the weight function 

(9.8) signx(x + x) v j(x 2 + a — c 2 ) 

on the union of two intervals [—\J ' /3 — a + c 2 , — c] U [c, \J '(3 — a + c 2 ] and satisfy the 
recurrence relation 

(9.9) S n+1 + (-l) n X Sn(x) + v n S n -!(x) = xS n (x), 
where v n are given by 

(9.10) V 2n = -C n , V 2n +1 = ~A n . 

The converse statement is also true. 

The orthogonality part of this Lemma can be checked by straightforward com- 
putations. Then, substituting (9.7) into the corresponding three-term recurrence 
relations we arrive at (9.9) and (9.10). For more details, see [13]. 

In what follows, we shall call the polynomials P n (x), the companion polynomials 
with respect to P n (x) (sometimes the polynomials P n {x) are referred to as the kernel 
polynomials [14]). 

10. A SCHEME FOR REDUCING THE POLYNOMIALS Q n (x; X) TO A-INDEPENDENT 

ORTHOGONAL POLYNOMIALS 

In this Section, we show how to reduce the polynomials Q n {x; A) defined by (7.3) 
to a A-independent system of orthogonal polynomials. 

First, let us apply Lemma 9.1 to Q n (x; A) for 9 = A — 1 and 9 = A + 1. Namely, 
if we put 9 = A — 1 then we can deduce that 
(10.1) 

-A(l + a„_i) if n is even 



-l + a n _i if n is odd 



Next, introducing the polynomials 

Q n +i(x; A) - A n Q n (x; A) 



^ _ ( — 1 — a n if n is even _ f 

" \ — A(l — a n ) if n is odd ' ™ \ 
Next, introducing the polynomi 

(10.2) Q n (x;X) 

x — A + 1 

we arrive at the the recurrence relation 

(10.3) Q n +i{x; A) + b n Q»(x; A) + u n Q n -i(x; A) = xQ n (x; A), 
where 

(10.4) &„ = (-l)"(A + l), u n = -A(l + a n )(l + (-l)"a n _i). 

Similarly, for 9 = A + 1 we have that 
(10.5) 

^ f 1 — a n if n is even f— A(l + a„_i) if n is even 

" \A(1 — a n ) if n is odd ' " \ — 1 — a n -i if n is odd 

Then the corresponding transformed polynomials Q„(x; A) have the recurrence co- 
efficients 

(10.6) &„ = (-l)"(A-l), u n = \(l + a n - 1 )(l-a n ). 
Clearly, both cases can be rewritten as follows 

(10.7) Q n +i(x; A) + (-l) n (diA + d )Q n (x; A) + A<g n _i(z; A) = xQ n (x; A), 
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where do, d\, and u* n > do not depend on A. 

The following Lemma reduces the polynomials Q n +\{x;X) to the polynomials 
obtained in Lemma 9.3. 

Lemma 10.1. Let the polynomials Q n (x;X) satisfy (10.7) and be orthogonal with 
respect to a weight function w\(x) on the set E\ for positive A. Then the monic 
polynomials defined via 

(10.8) S n (x; A) = (VX)- n Q n (V\x; A) 

are orthogonal with respect to the weight function w\(VXx) on the set E* x = {x : 
VXx ^ E\} and satisfy 

(10.9) S n+1 (x; A) + (-l) n X S n (x; A) + <S„_i(z; A) = xS n (x; A), 

where x — diy/X + d a . 
The converse is also true. 

The proof is straightforward by making the substitution x — > VXx. 

Now, we can apply Lemma 9.3 to represent the polynomials S n (x; A) in terms of 
a A-indcpcndcnt system of orthogonal polynomials. Also, Lemma 9.3 gives us the 
explicit formula for the weight function of the polynomials S n (x; A). 

Summing up, let us notice that the scheme given in this Section gives us the 
possibility to reduce the A-dependent system of polynomials Q n (x; A) to a system 
of A-indcpendent orthogonal polynomials P n (x) given by (9.1). Namely, we reduced 
Q n (x\ A) to the polynomials Q n (x; A) satisfying (10.7). Then Lemma 10.1 leads us 
to the polynomials S n , which according to Lemma 9.3 can be constructed from the 
polynomials P n (x) given by (9.1). 

Finally, it is worth noting that Lemma 9.3 and Lemma 10.1 give an efficient way 
to construct explicit families of Jacobi matrices of the form 

Ji + AJ 2 , 

where J\ and Ji are bidiagonal matrices. The corresponding eigenvalue problems 
with spectral parameter x lead to orthogonal polynomials and to Laurent biorthog- 
onal polynomials if we fix x and consider the problem with respect to A [26], [41]. 

11. From the Jacobi OPUC to the big -1 Jacobi polynomials 

To fit the big -1 Jacobi polynomials in this context, let us consider the polyno- 
mials 

(11.1) Q* n (x;X)=g- n Q n (gx;X). 
Then (7.3) yields the following 

(11.2) Q* n+ i(x; X) + bl(X)Q* n (x- X) + <(A)Q;_ 1 (x; A) - xQ* n (x; A), 
where the coefficients are as follows 

(11.3) b* = ±(\-l-A n -C n ), u* n = ^C n A n . 1 , 
and A n , C n are defined in (10.1). 
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Next we compute the coefficients A n and C n for the reflection parameters a n 
given by (6.2). We have 

(11.4) ' 

A = J if » is ™ f if „ is even 

" if » is odd ' " ^ » * odd' 

Set A = and 5 = — Now, introducing the new parameters 

(11.5) A' n = ±A n , C' n = -C n 
and choosing the parametrization 

(11-6) ^ = — ' 77 = —' 

we can easily see that the coefficients A' n and C' n correspond to those of the big 
-1 Jacobi polynomials introduced in [38] with a shifted spectral parameter. Recall 
that the big -1 Jacobi polynomials are polynomials orthogonal on the union of two 
intervals [— 1, — c] U [c, 1] with respect to the weight function 

(11.7) w^- x \x) = signx(x + l)(x - c)(l - x 2 )^-^' 2 {x 2 - c 2 )^ 1 )/ 2 . 

Notice that the case A = 1 corresponds to c = 0. This special case corresponds to 
the little -1 Jacobi polynomials. In this instance the two intervals of orthogonality 
connect to a form the single interval [—1,1]. 

The big -1 Jacobi polynomials P^ a '^\x; c) corresponding to the weight function 

(11.7) are the eigensolutions of [38] 

(11.8) LP^\x;c)=X n P^(x;c), 
where the operator 

(11.9) - g Q (x)(R -I)+ 9l {x)d x R 
with 

(a + /3 + l)x 2 + (ca - /3)x + c 2{x - \){x + c) 

(11.10) g {x) = 5 , gi(x) = , 

x z x 

I is the identity operator and R the reflection operator Rf(x) = f(—x). 
The eigenvalues of L^ a ^^ are 

2n, n even 
-2(a + P + n+ 1), n odd ' 

As shown in [39] the operator (11.9) is the most general operator of the first order 
in d x and involving R that has orthogonal polynomials as eigenfunctions. 

We thus see that the orthogonal polynomials Q n (x; A) obtained by the SDG map 
from the Jacobi OPUC, satisfy a Dunkl type eigenvalue equation. This equation 
can easily be obtained from (11.8) by an affine transformation of the argument 
x — > j-ix + v with appropriately chosen constants ji, v. 

Note that when £ = t] = —1/2 in (6.2), a n = 0, n = 0, 1, 2, . . . , i.e. this case is 
equivalent to the 1-periodic case considered in Section 8. This corresponds to the 
case a = j3 = of the big -1 Jacobi polynomials. In particular, the Chebyshev 
polynomials of the third kind [29] (corresponding to A = 1, or equivalently, c = 0) 

cos(n(r + l/2)) 

V n (x) = -— , X = COST 

cos(r/2) 



(11.11) A„ = 
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coincide with the little -1 Jacobi polynomials [37] with a = j3 = and hence they 
satisfy the Dunkl type eigenvalue equation 

(11.12) 2(x - l)V;(-x) + V n (-x) = (-l) n (2n + l)V n (x). 

Similarly, for the Chebyshev polynomials of the fourth kind [29] 

sin(n(r + 1/2)) 

W n (x) = . -. ... , x = cosr 

sin(r/2) 

one has the eigenvalue equation 

(11.13) 2{x + l)W' n (-x) + W n {-x) = (-l) n (2n+l)W„(a;). 

In fact, equation (11.13) is a simple consequence of (11.12) because of the property 
V n (-x) = (-l) n W n (x) [29]. 
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